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ABSTRACT 



Context. The deuterium fractionation, Dfrac, has been proposed as an evolutionary indicator in pre-protostellar and protostellar cores 
of low-mass star-forming regions. 

Aims. We investigate Df^ac, with high angular resolution, in the cluster environment surrounding the UCH 11 region IRAS 20293+3952. 
Methods. We performed high angular resolution observations with the IRAM Plateau de Bure Interferometer (PdBI) of the ortho- 
NH2D lii-loi line at 85.926 GHz and compared them with previously reported VLA NH3 data. 

Results. We detected strong NH2D emission toward the pre-protostellar cores identified in NH3 and dust emission, all located in 
the vicinity of the UCH 11 region IRAS 20293+3952. We found high values of Dfrac -0.1-0.8 in all the pre-protostellar cores and low 
values, i)frac< 0.1, associated with young stellar objects. 

Conclusions. The high values of Df^ac in pre-protostellar cores could be indicative of evolution, although outflow interactions and UV 
radiation could also play a role. 
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1. Introduction 

Characterizing the diflTerent evolutionary stages before and after 
] the formation of a star is crucial for understanding the process 

■ of star formation itself. Dense cores, where stars are born, are 
[ mainly studied through molecular emission of dense gas tracers, 

and N-bearing molecules are widely used because they do not 
freeze out on to dust grains until very high densities are reached 
(~ 10^ cm~^: [Flower et al.ll2006h . Several column density ratios 
. of molecules tracing dense cores have been prop osed as good 

■ cheniical clocks, such as NH3/N2H+ or CN/N2H+ (iHotzel et alJ 
; l2004l:lFuente et al.ll2005h . Among these, the ratio of a deuterated 

■ species over its counterpart containing H, i. e., the deuterium 

■ fractionation Dfrac, has been found to be a good tracer of the 
evolutionary stage of dense cores. Both observations and mod- 
els strongly suggest that in cold (7 < 20 K) and dense cores (n ^ 
10^ cm"^) C-bearing molecules are expected to deplete onto dust 
grain s, leading to an enhancement of the deuterium fractionation 
fe^ [Roberts & Millaill2Q0Ql: iBacmann et al.ll2003l: iPillai et al.1 
1200 7*) because the H2D^ ion, the progenitor of most of the 
deuterated species, including NH2D, is not destroyed by CO. 

In particular, in low-mass star-forming regions, Df^ac is found 
to increas e until the onset of star formation and to decr ease af- 
terwards (ICrapsi et al.ll2005l: lEmprechtinger et al.ll2009l) . In the 
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high-mass regime, Chen e t al. 1 120T0I) performed sinsle-dish ob- 
servations toward massive dense cores and find a decreasing 
trend of Df^ac with the evolutionary stage during the protostel- 
lar phase. However, since massive star formation takes place 
in cluster environments, to study Z)frac one needs to carry out 
observations with high angular resolution in order to separate 
the emission from each individual core. So far, interferomet- 
ric observations of deuterated specie s in massive star-forming 
regions have only been reported by iFontani eta l. (2008) and 
[Sandell & WrighH (2010), but only the work of Fontani et all 
(l2QQ8h evaluates Z)frac in the region, where they find two dense 
cores with Z)frac- O-l- 

In this letter we present high angular resolution observations 
of the ortho-NH2D In-loi line at 85.926 GHz carried out with 
the Plateau de Bure Interferometer toward the high-mass star- 
forming region IRAS 20293+3952, located at 2 kpc of distance 
(lBeutheret al.ll2004bl) and with a luminosity of 6300 Lq. The 
region is associated with an UCH 11 region at the border of a 
cloud of dense ga s mapped in NH3 with high angular resolution 
(iPalau et al.ll2007l) . and it harbors a rich variety of young stellar 
objects (YSOs) and dense cores at diflferent evolutionary stages, 
being thus an excellent target to study Z)frac- While the north- 
ern side of the main cloud of the region c ontains several YSOs 
driving at least f our molecular outflows (iBeuther et al.ll2004al: 
iPalau et all l2007l) . the southern side is mainly populated with 
starless core candidates (BIMA 3 and BIMA 4, two faint mil- 
limeter condensations). In addition, there is a small cloud to the 
northwest, the western cloud, which seems to harbor a very com- 
pact starless core (see Fig.[T]for an overview of the region). The 
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Fig. 1. Contours: NH2D zero-order moment. Contours start 
at 6 %, increasing in steps of 10 % of the peak intensity, 
0.58 J y beam~^ k m s"^ Color scale: NH3 (1,1) zero-order mo- 
ment (iPalau et al.ll2007l) . The NH2D synthesized beam, 4732 x 
3705 (P. A.= 111°), and the NH3 synthesized beam are shown 
in the bottom left and bottom right corners, respectively. White 
star: position of the UCHii region (IRAS 20293-h3952). Plus 
signs: YSOs embedded in the main cloud. Black triangles: NH3 
column density peaks in the starless cores BIMA 3, BIMA 4, and 
the western cloud. Main cores of NH2D are labeled as western 
cloud, dust ridge, BIMA 2-S, and BIMA 4. The blue and red con- 
tinuum/dashed arrows mark the direction of powerful/less pow- 
erful molecular outflows, respectiv ely (Beuther et al. 2004a ). 
Outflow B is deflected by BIMA 4 (iPalau et al.ll2007l) and the 
dashed lines indicate the high- velocity gas around BIMA 4. 



current NH2D observations allowed us to estimate for the first 
time Dfrac from NH2D/NH3with high angular resolution toward 
a massive star-forming region containing dense cores and YSOs 
at difl'erent evolutionary stages. 

2, Observations 

The IRAM Plateau de Bure Interferometer (PdBI) was used to 
observe the NH2D In - loi molecular transition at 85.926 GHz 
toward IRAS 20293-1-3952. The observations were carried out 
during 2008 June 16 and December 4, with the array in the 
D (4 antennas) and C (6 antennas) configurations, respectively, 
providing projected baselines ranging from 17.5 m to 175 m. 
The phase center was a(J2000) = 20^31^12^.7, ^(J2000) = 
-f40°03'1374. The typical system temperatures for the receivers 
at 3 mm were -150-300 K during June 16, and -100-150 K 
during December 4. Atmospheric phase correction was applied. 

The 3 mm receiver was tuned at 86.75433 GHz in the lower 
sideband. The spectral setup included other molecular transi- 
tions, which will be presented in a subsequent paper (Busquet 
et al., in prep). For the ortho-NH2D, we used a correlator unit 
of 20 MHz of bandwidth and 512 spectral channels, which pro- 
vides a spectral resolution of ~ 0.039 MHz (~ 0.14 km s"^). 
The FWHM of the primary beam at 3 mm was ~ 56''. Bandpass 
calibration was performed by observing the quasars 3C454.3 on 



June 16 and 3C273 on December 4. Amplitude and phase cali- 
brations were achieved by monitoring MWC 349 and 2005-1-403 
for both days. The phase rms was ~ 20°. The absolute flux den- 
sity scale was determined from MWC 349, with an estimated 
uncertainty ~ 10%. The data were calibrated with the program 
CLIC and imaged with MAPPING, and both are part of the 
GILDAS software package. Imaging was performed with natu- 
ral weighting, obtaining a synthesized beam of 4V32 x 3V05 with 
P. A.= 111°, and rms of 20.7 mJy beam"^ per channel. 

3. Results 

In Figd] we present the NH2D In-l pi zero-ord er map integrated 
for the 6 hyperfine transitions (see lOlberg et al. 1985 for a de- 
scrip tion of the transitio n) overlaid onto the NH3 (1,1) emission 
from lPalau et al.l (l2007h . The integrated NH2D emission toward 
the western cloud presents a compact morphology with a de- 
convolved size of 6r5 x 376 (~ 0.06 pc). We also found NH2D 
emission associated with the southern side of the main cloud, 
covering a spatial extension in the east- west direction of ~ 33'' 
(0.33 pc). The emission consists of three main cores, one clearly 
associated with BIMA 4, another NH2D core coincident with the 
dust ridge, and a third core located ~ 7" south of BIMA 2 (here- 
after BIMA 2-S). In addition, we detected faint NH2D emission 
associated with BIMA 3. It is worth noting that NH2D is not (or 
marginally) detected on the northern side of the main cloud. In 
this letter we follo w the nomenclatur e for the YSOs and starless 
core candidates of Palau et al.l (|2007'). 

We studied the deuterium fractionation of the region by 
first computing the column density map of NH2D and then the 
A/^(NH2D)/A/^(NH3) column density ratio. For this, we extracted 
the spectra for positions on a grid of 2'' x2", a nd adopted the hy- 
perfine frequencies listed in llTine et al.l (l2000l) . The column den- 
sities are only reported for spectra with a peak intensity greater 
than 5 cr, and with all the hyperfine components detected. 

Assuming that all levels are populated according to the same 
excitation temperature, Tex, the column density of the asymmet- 
ric top molecule NH2D is given by 

[=^1 — "■[oferiM'^-dii^n 



where Av is the linewidth, the optical depth of the main 
line (lii-loi, F=2-2), v the frequency of the transition (y = 
85.926 GHz), Eu the energy of the upper level (Eu = 20.68 K), 
Aui = 5.8637 X 10"^ s"^ is the Einstein coefficient, and 
/^(T^ex) = (/zy/^)/(^^^/^^- - 1), where h and k are the Planck and 
Boltzmann constants, respectively . From the Cologne Database 
for Molecular Spectroscopy (CDMS; Miiller et al. 2001), the up- 
per level degeneracy is g„ = 15, and the partition function Q(Tq^) 
is estimated as Q(Tq^)= a -h I3TqJ^^, where a = 3.899 and 
j3 = 0.751 are the best-fit parameters from a fit to the partition 
function at the diff'erent temperatures giv en in CDMS. The stan- 
dard ortho/para ratio of 3 (e. g., Walmsley et al.lll98"7l : lTine et"aD 
I2OOO ) is already included in the partition function. 

Figure [2] shows a map of the linewidth (Av) derived from the 
hyperfine fit and the total column density A/^(NH2D) toward the 
southern side of the main cloud (i. e., dust ridge core, BIMA 2-S, 
BIMA 3, and BIMA 4, hereafter referred to as the main cloud; 
see left panels) and the western cloud (right panels). The prop- 
erties of the western cloud and the main cloud are somewhat 
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Fig. 2. Color-scale maps of linewidth and A/^(NH2D) overlaid on the NH2D 1 n-loi integrated emission (contours) for the main cloud 
(left panels) and the western cloud (right panels). Triangles mark the NH3 column density peaks of BIMA 3, BIMA 4 (left), and the 
western cloud (right). In all panels the synthesized beam is shown in the bottom right corner. 
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Fig.3. Spectra (black line) toward 5 positions of the IRAS 20293+3952 for NH2D In-loi (top) and NH3 (1,1) (bottom). The 5 
positions, which are labeled in the top right corner of each panel, are, from left to right WC (western cloud), BIMA 2-S, DR (dust 
ridge core), BIMA 3, and BIMA 4. The red lines show the fit to the hyperfine structure obtained as illustrated in Sect. 3. 



Table 1. Summary of the main physical properties of selected cores in IRAS 20293+3952 





Qf(J2000) 


(5(J2000) 


^ ex 




7- c 
^ rot 


A^(NH2D)^ 


//(NHg)^ 




NH3/ 


Evol. 


Core 


(h m s) 


C ' ") 


(K) 


' m 


(K) 


(xlO^^ cm-2) 


(xlO^^ cm-2) 


frac 




stage^ 


Western cloud 


20:31:10.70 


40:03:28.2 


4.8+0.2 


1.5 


15+3 


10+ 2 


14+ 2 


0.8+0.2 


300 


PPC 


BIMA 2-S 


20:31:12.34 


40:03:05.7 


5.4+0.3 


1.1 


19+3 


2.8+0.6 


6.7+0.8 


0.4+0.1 


90 


PPC 


Dust ridge 


20:31:12.98 


40:03:08.5 


4.6+0.3 


1.0 


16+3 


2.8+0.7 


4.9+0.9 


0.6+0.2 


90 


PPC 


BIMA 3 


20:31:13.92 


40:03:00.9 


5.4+1.0 


0.4 


14+2 


1.9+0.5 


19+ 3 


0.10+0.04 


300 


PPC 


BIMA 4 


20:31:14.56 


40:03:06.8 


7.5+1.6 


0.3 


14+2 


1.2+0.4 


7.4+0.8 


0.16+0.05 


200 


PPC 


BIMAl 


20:31:12.77 


40:03:22.6 






24+3 


< LIS 


19+3 


<0.06 


50 


0, no IR 


IRS 5 


20:31:13.41 


40:03:13.7 






24+3 


< 1.52 


11+2 


<0.1 


50 


IR 



Notes. Tex derived from the ou tput parameters of the NH2D ln-loi hyperfine fit. Derived from the fits to the hyperfine structure of 
NH2D lii-loi. Obtained from P alau etal.1 (2007). Column densities derived using the same uv range and convolving the images to a 
circular beam of 1" . D^slc = A/^(NH2D)/A/^(NH3). PPC: pre-protostellar core; O: molecular outflow. 3cr upper limit estimated adopting 
A V from NH3 and T^^= 5 K. 



different. The total optical depth (ti^^-Iq^ = 2t^) in the west- 
ern cloud is in the range ~ 4-15, with a typical uncertainty of 
-1.1, whereas the total optical depth is lower toward the main 
cloud, between 0.5 and 4, with a typical uncertainty of ~ 0.5. 
Concerning the linewidth in these clouds (see Fig. Oa), we find 
line broadening (Av^ 0.75-0.9 km s"^) associated with the west- 
ern cloud. Toward BIMA 3 and BIMA 4, we also find high val- 
ues for the linewidth, around ~ 0.8-1.2 km s"\ while the dust 



ridge core and BIMA 2-S appear to be more quiescent, with 
linewidths of -0.5-0.7 km s"^. Since it has been suggested that 
one of the outflows of the region, outflow B, is interacting with 
BIMA 4 (Palau et al. 2007), this could produce the fine broad- 
ening seen in NH2D in BIMA 3 and BIMA 4. Additionally, the 
line broadening found in the western cloud is spatially coinci- 
dent ^in_2roiection}w a high- velocity feature of outflow D 
(Beuther et al.l l2QQ4ab . Thus, the broad lines seen toward the 
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western cloud, BIMA 3, and BIMA 4, suggest that the deuter- 
ated gas is being perturbed by the passage of the outflows (see 
Fig.©. 

The NH2D column density, corrected for the primary beam 
response, also presents significant diff'erences between the west- 
ern cloud and the main cloud. While the column density is, on 
average, A/^(NH2D)^ 25 x 10^^ cm"^ in the western cloud, in the 
main cloud the column density is slightly lower, in the range 1- 
8x10^"^ cm"^, reaching its maximum value of 8 x 10^"^ cm"^ close 
to the peak position of BIMA 4 (see Fig.|2]b). The uncertainty of 
the NH2D column density is -25-35 %. 

In order to properly estimate the deuterium fractionation, de- 
fined as Dfrac = A^(NH2D)/A^(NH3), we made the NH3 VLA im- 
ages using the same uv range as the PdBI data (5-50 kA), estimat- 
ing the column densities for both NH3 and NH2D for the same 
angular scales. Finally, we convolved the NH2D and NH3 emis- 
sion to a circular beam of 7'' (the major axis of the NH3 beam). 
In Fig. [3] we present the spectra obtained at the NH2D emission 
peak of each condensation (i. e., western cloud, BIMA 2-S, dust 
ridge core, BIMA 3, and BIMA 4) together with the hyperfine fit 
obtained toward these positions. In Table[T]we list the excitation 
temperature, Tex, the rotational temperature, T^ou the NH2D and 
NH3 column densities, Z)frac, and the NH3/N2H^ ratio for each 
core and a few YSOs. Toward the YSOs BIMA 1 and IRS 5, 
we report on upper limits, with Z)frac< 0-1, and we cannot draw 
any conclusion for the behavior of Df^^c in the protostellar phase. 
More interestingly, in the western cloud Df^ac is ~ 0.8, which is 
the highest value of Df^ac in the region and among the highest re- 
ported in the literature (e. g. JCrapsi et al .120071: IPiUai et al.l2007t 
iFontani etar.11200 8). In the main cloud, Dfrac presents significant 
variations among the diff'erent cores, with Df^ac decreasing from 
the northwest (Z)frac-0.5 in BIMA 2-S and the dust ridge core) to 
the southeast (Dfrac-O.l in BIMA 3 and BIMA 4). This suggests 
a chemical diff'erentiation along the main cloud. 

4. Discussion and summary 

Our high angular resolution study of the NH2D toward the 
massive star-forming region IRAS 20293 -h3952 reveals strong 
NH2D emission toward starless cores, whereas NH2D is not 
(or marginally) detected in cores containing YSOs, which sug- 
gests that the production of NH2D is more eff^ecti ve in the 
pre-pr otostellar phase than in the protostellar phase. Pata u et aP 
(I2007h notice that the starless cores in this region seem to be 
predominant on the southern side of the main cloud and in the 
western cloud, while the northern side of the main cloud har- 
bors all the YSOs known in the region, suggestive of the dense 
gas in the main cloud being progressively more evolved as it 
moves from south to north. In addition, chemical diff'erentiation 
among pre-protost ellar and protostellar cores was also found by 
iPalau et al.l (12007 *) using the NH3/N2H+ ratio, which was high 
for pre-protstellar cores and low in protostellar cores (see col. 9 
of Table[T]). Thus, for this region, the behavior of Z)frac, measured 
from NH2D/NH3, is similar to the behavior of NH3/N2H^ ratio, 
suggesting that both ratios can be used to distinguish between 
pre-protostellar and protostellar cores and that both ratios could 
be related with the evolutionary stage of the dense gas. 

A possible interpretation of the diff'erences in Df^ac seen 
in the pre-protostellar cores of region could be that they are 
in diff'erent evolutionary stages. According to the study of 
ICrapsi et aP (|2005). there is an increasing trend for Dfrac as 
the starless core approaches the onset of gravitational collapse 
(from 0.03-0.1 in the youngest cores to 0.1-0.4 toward the most 
evolved cores). This would indicate that the western cloud is the 



most evolved pre-protostellar core and that BIMA 3 and BIMA 4 
are less evolved. However, in regions of massive star formation, 
typically associated with clustered environments, other factors, 
like temperature, UV radiation, and/or molecular outflows, can 
play important roles in altering the chemistry, and then it is not 
clear whether this trend is related to the evolutionary stage of 
pre-protostellar cores. 

A comparison of Df^ac between the western cloud and 
BIMA 3/BIMA 4, all of them having similar temperatures but 
very difl'erent values of Dfrac (see Table [TJ, indicates that tem- 
perature is not an imp ortant factor in dete rmining Df^ac for tem- 
peratures around 15 K. iPalau et al.l (12007) find evidence that UV 
radiation from the UCHii region aff'ects the chemistry at the 
western edge of the main cloud facing the UCH 11 region, so it 
could aff'ect BIMA 2-S and the western cloud. In particular, the 
presence of a cavity between the UCH 11 region and the western 
cloud suggests that this cloud could be photo-illuminated by the 
UCH II region. However, while the western cloud and BIMA 2- 
S could be locally aff'ected by the UV radiation, this is not the 
case for BIMA 3 and BIMA 4, for which the high visual extinc- 
tion in the main cloud prevents UV photons from penetrating. 
Finally, Font ani et al.l (2009) print out that shocks in outflows 
could modify Z)frac- This could be the case for the western cloud, 
BIMA 3, and for BIMA 4. While the interaction of an outflow 
with the western cloud is not evident, Patau et at.i (|2007) have 
already proposed that the powerful outflow B is interacting with 
BIMA 4, producing the deflection of the outflow and thus the 
ejection of high- velocity material in different directions. This 
could afl'ect BIMA 4 and BIMA 3 beca use high- velocity SiO 
is seen close to these cores (see Fig.fTland lBeuther et al.ll2004ah . 
This interpretation is reinforced by the line broadening of NH2D 
observed towards BIMA 3 and BIMA 4 (see Fig.[2|. In conclu- 
sion, while Dfrac in IRAS 20293 -h 3 95 2 may be locally afl'ected 
by the interaction of outflows and UV radiation, Df^ac is lower 
toward YSOs than toward pre-protostellar cores, with a possible 
evolutionary trend in the pre-protostellar phase, which deserves 
further study. 
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